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Abstract. Fluctuation and correlation observables are often measured using multi-particle 
correlation methods and therefore mutually probe the origins of genuine correlations present in 
multi-particle distribution functions. We investigate the common influence of correlations arising 
from the spatially inhomogeneous initial state on multiplicity and momentum fluctuations as 
well as flow fluctuations. Although these observables reflect different aspects of the initial 
state, taken together, they can constrain a correlation scale set at the earliest moments of the 
collision. We calculate both the correlation scale in an initial stage Glasma flux tube picture and 
the modification to these correlations from later stage hydrodynamic fiow and find quantitative 
agreement with experimental measurements over a range of collision systems and energies. 



1. Introduction 

Fluctuating initial conditions of relativistic heavy ion collisions influence not only flow and but 
also fluctuation measurements. The random distribution of transverse parton density "lumps" 
results in arbitrary event shapes. Harmonic flow is generated largely by the global spatial 
anisotropy of the initial shape and induces a global correlation in the final momenta of measured 
hadrons. Conversely, the mere existence of parton density lumps results in local correlations; 
partons emerging from the same source are more likely to share the same phase space throughout 
the collision lifetime. We argue that the same local correlations that give rise to multiplicity 
and pt fluctuations drive fluctuations of the harmonic flow coefficients [U |2j . 

To begin, in Sec. [2]we describe multiplicity, pt, and Vn fluctuations as well as the ridge in terms 
of two-particle correlations. Using a cumulant expansion of the pair distribution reminiscent of 
those in Refs. [3 1^, we show that all of the observables are linked by a common correlation 
function. The normalization of this correlation function sets a common scale that governs most 
of the centrality and center of mass collision energy dependence. 

In Sec. [3] we investigate contributions to the correlation function from transversely local 
longitudinally long range correlations. We build on an approach started in Refs. [5l |6] in 
which long range correlations result from the fragmentation of Glasma flux tubes. In this 
formulation local spatial correlations emerge from fluctuations in the number and distribution 
of flux tubes. These spatial correlations are then modified by transverse expansion, giving rise 
to azimuthal correlations. A similar physical picture motivates studies using a wide range of 
different techniques [TIElEllinKIIlOElE]. 



Finally in Sec. |4]we discuss how fluctuation measurements calculated in a unified framework 
can be used to constrain theories of the initial state of heavy ion collisions. 



2. Correlations, Fluctuations, and Flow 

Correlation measurements largely rely on the distribution of hadron pairs P2(pijP2) = 
(iA^/(ipi(ip2. Using the cumulant expansion method the pair distribution can be written as 

P2(Pl,P2) = Pl(Pl)Pl(P2) +r(pi,P2), (1) 

where pi{p) = dN/dp is the singles distribution, and r(pi,p2) represents the genuine two- 
particle correlation function [U O El Uj . 

In this section we outline the role of correlations, T(pi,p2), in determining fluctuation and 
flow observables and leave our model description to Sec. [3j In the absence of correlations, the 
pair distribution factorizes into the square of the singles distribution, /02(Pi; P2) Pi(Pi)pi(P2)- 
This assumption is often cited as a signature of harmonic flow, given that the event geometry 
is the dominant contributor to Vn measurements. It is informative, however, to discuss the 
implications of relaxing this assumption when considering flow fluctuations. To address this 
we study other fluctuation observables such as multiplicity and momentum fluctuations which 
do not rely on event geometry but still reflect the existence of genuine correlations. We then 
find that although harmonic flow coefficients do depend on the global correlations induced by 
geometry, flow fluctuations emerge largely from the the same correlations as other fluctuation 
observables. 



2.1. Multiplicity Fluctuations 

The existence of genuine correlations in the pair distribution signifies non-Poissonian particle 
production; all particles are not emitted from random independent sources Rearranging 
as 

'^(Pl, P2) = /02(Pl, P2) - Pl(Pl)Pl(P2), (2) 

one can define the " robust variance" , 

1 /■ . ^ , , (N(N- 1)) - (iV)2 
^=^J r{pi,P2)dpidp2 = ^J^, (3) 

We can rewrite the right most term of Q as 7^ = Var{N) — {N)/{N)'^. Possonian particle 
production requires Var{N) = {N), and consequently, non-zero values of TZ indicate a correlation 
in the production of hadrons. We will therefore refer to TZ as the correlation strength in this 
work. Significantly, the PHENIX collaboration measures multiplicity fluctuations using the 
width of the negative binomial distribution of produced hadrons |16j . The inverse of the width 
parameter is equivalent to the correlation strength, TZ = kj^^^ [T7]. The results are shown in 
the bottom panel of Fig. [T] This quantity sets the scale of all fluctuation measurements using 
two-particle correlations. 

In Sec. [3] we will argue that these correlations originate from the fluctuating initial parton 
density. Pairs produced together in the hot spot or density lump are more likely to be correlated. 
We stress here that the shape of the event due to the distribution of hot spots does not contribute 
to multiplicity fluctuations. Integration of ^ over all momenta, including azimuthal angles, 
sums all correlated particle pairs, regardless of direction. 



2.2. Transverse Momentum Fluctuations 

Comparably to multiplicity fluctuations, many describe the dynamic fluctuations of transverse 
momentum using the covariance 



{SmSp,2) = . (4) 

where 6pti = pu — {pt) and the average transverse momentum is {pt) = (Pt) / {N) for Pt = J2iPti 
the total momentum in an event \18\ [191 [20l [2T] . This quantity vanishes when particles i and j 
are uncorrelated. Equation Q sums pairs of differences in particle momenta from the average 
and can be written in terms of the pair correlation function ^ as 

{SptiSpt2) = J dpidp2 J^j^^~]yj^Pi'^^Pi'^- (^) 

Pairs in which both particles have higher than average momentum add to {dptiSpt2)- Lower- 
than-average pairs also add to the covariance, while high/low pairs subtract from it. In global 
equilibrium {6ptiSpt2) = 0. The presence of hot spots makes {6pti6pt2) > (as would cold spots) 
|22j . Motion of the sources would further enhance this quantity It follows that both jets 
and flow add to the pt covariance. 

It is important to notice that momentum fluctuations depend on the same correlations as 
multiplicity fluctuations, the only difference being momentum weights. Like ([S]) the integration in 
([5| counts pairs independent of their direction and therefore does not depend on the asymmetry 
of the system. Because of the momentum weights, ([S]) is sensitive to the average expansion and 
temperature and reflects the degree of thermalization of the system. 

2.3. Flow Fluctuations 

Unlike multiplicity and momentum fluctuations, harmonic flow correlation and fluctuation 
measurements depend on the azimuthal anisotropy of the collision system. The anisotropy 
of the collision can be characterized by the moments 

(v„) = (cosn((/> - Vbp)>> (6) 

where (...) represents an average over particles and events and ip^p is the reaction plane 
angle defined by the plane spanned by the impact parameter b and the beam direction. 
Since the reaction plane cannot be measured, experiments turn to the two-particle correlation 
measurement 

r .2 _ (E.^jCOsn((/>, -</.,■)) 
"""^^^ " (A^(iV-l)) ^ ' 

[31 in [231 1211 [25] ■ Being a two-particle correlation, we can cast ([7]) in terms of the pair distribution 

vn{2}' = J ^0|^^cosn(<^i -<A2)cipidp2. (8) 



Using ([T]), ^ can be written as 

t;„{2}2 = {v^f + 2al (9) 

where the flow coefficient relative to the reaction plane is (u„) = / pi(p) cosn(i;^ — ^^p)(ip 
comes from the single particle distribution and the factor of two in the fluctuation term a is 
conventional [26j. The effects of geometry primarily enter from the anisotropy of the single 



particle distribution characterized by ([6]). Genuine correlations contributing to ([S]) and ([S]) 
contribute to flow fluctuations 



(Tn = J dpidp2 2(jv(Ar- i)) ^'^^ 

Ulike ^ and ([5]), ([lO]) sums the relative azimuthal angles A(j) = 0i — 02 of correlated pairs. 
Anisotropic flow would modify correlated pairs emerging from common density sources since 
those pairs would experience a different push depending on their location in the transverse 
plane. 

Thus far we have only mentioned correlations among flowing particles. The correlation 
function r(pi,p2) also contains correlations emerging from jets, resonance decays, HBT effect, 
etc.. These correlations may or may not be modified by flow but contribute to all mentioned 
observables. To eliminate these effects we can consider long range correlations and measurements 
of observables with rapidity separations. 

2.4. The Ridge 

Ridge measurements study two-particle azimuthal correlations at small and large 
(pseudo)rapidity separations Arj. Correlations at rapidity separations larger than ~ 1 — 2 
units cannot come from jets, resonance decays or other short range phenomenon and causality 
dictates that such long range correlations originate in the early stages of the collision |271 [5], 

The soft ridge observable developed by STAR can be shown to relate directly to flow 
coefficients as well as fiow fiuctuations by 



= 0:= ^ 2 E (^") + : \ ■ (11) 



[2]. Here the (vn) terms follow from ([g]), and pref is the experimental mixed-event background. 
We can isolate the contribution of flow fluctuations using 



dN 2 f Ap{A(l),Arj) 



cos(nA(/)) dA4>dAr], (12) 

FS 



where {Ap{A(j), Ar])/ ^p,.c{)\fs is the flow-subtracted ridge amplitude on the near side, centered 
at A(/) = At/ = as reported by STAR [2811291 [30]. 

Since ridge measurements rely on both flow and fluctuations, they can be compared to 
independent flow and flow fluctuation measurements to determine the relative contributions 
of long and short range correlations. Moreover, calculating all of the observables mentioned in 
Sec. [2] within the same framework can constrain the common correlation function ^ present in 
all measurements. 



3. Glasma Correlations and Results 

We argue that multiplicity, momentum, and flow fluctuations originate from initial state spatial 
correlations. These correlations emerge from fluctuating parton density distributions. Density 
lumps or hot spots emerge reflecting the positions of partons sources such as wounded nucleons 
or, in our case. Color Glass Condensate - Glasma flux tubes. We take correlated particles as 
those originating from the same source. 
The spatial correlation function 



c(xi,X2) = n2(xi,X2) - ni(xi)ni(x2) 



(13) 
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Figure 1. Prediction for IZdN/dy 
function on the number of 
participants A'part at three beam 
energies (top). Calculated TZ 
compared to PHENIX data from 
[16J (bottom). 



identifies the initial distribution of correlated pairs in the transverse space. If the pair 
distribution, n2(xi,X2), contains no correlations, it factorizes into the square of the singles 



distribution, ni(xj) and (13) vanishes. We write (13) in a Glasma flux tube picture as 



c(xi,X2) = (iV)27e5(rt)p^^(Rt). (14) 

The delta function in relative transverse elliptical coordinates = rti— rf2 enforces the condition 
that correlated partons originate from the same flux tube presuming the flux tube transverse size 
is small compared to that of the collision area. The flux tube probability distribution, written 
in average coordinates = (r^i + r(2)/2, is 

P.T(Rt)«^fl-^|- (15) 



Here the shape of (15) resembles the nuclear thickness function and vri?j^ is the area of the 



overlap region. The distribution (15) represents an average over all possible shapes, which we 
have taken to be a simple ellipse. 

To compute the correlation strength 7^, we imagine each event produces K flux tubes with 
transverse size ~ Q^^- saturation regime K is proportional to the transverse area 

divided by the area per flux tube, j31j . Allowing K to fluctuate from event to event with 
average {K), we calculate in Ref. [5] that 



(iV2) - (iV)2 _ (AT) 



Each Glasma flux tube yields an average multiplicity of ~ ^{Qs) gluons and as in Ref. [31] , 
the number of gluons in a rapidity interval Ay is then 

(iV) = {dN/dy)Ay ~ as-\Qs){K). (17) 

Finally, the Glasma correlation scale, 

TZdN /dy = na~^{Ql), (18) 

follows from the multiplication of ([s]) and (17). Notice ( |18[ ) is dimensionless and depends only 
on the saturation scale, which can be calculated from first principles. Measurements of the 
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Figure 2. Transverse momentum 
fluctuations {5pti5pt2) dN / drj 
function on the number of partici- 
pants A'^part at three beam energies. 
Data is from [211132] . 
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Figure 3. Measured f2{2} and 
t>2{4} from STAR [33] compared 
to calculations. The difference 
between the calculated curves is 
due to flow fluctuations. 



ridge at various beam energies, target masses, and centralities fix the dimensionless coefficient 
K and are in excellent accord with the leading-order dependence [51 [B] . 

Intrinsically correlated partons originate from the same transverse position and experience, 
on average, the same transverse momentum modulation from fiow. We represent the genuine 
two-particle momentum space correlation as 



r{pi-,P2) =j c(xi,X2)/(xi,pi)/(x2,p2)(iridr2 



(19) 



We evaluate (19) using a blast- wave description of the Boltzmann phase-space density /(x, p) 
on a Cooper- Frye freeze out surface dP [Il[2]. Notice that full integration of (19) yields {N)'^TZ 



as does (13), tying azimuthal correlations to multiplicity and transverse momentum fiuctuations 



Using (19) we calculate ([S]) and find the results shown in Fig. [2j We find good agreement 
in central collisions where the assumption of gluon saturation necessary for the production of 
Glasma flux tubes is most valid. Moreover, the onset of thermalization in peripheral collisions 
should modify the correlation function. In particular, this effect has been shown to modify 
{6ptiSpt2) at low numbers of participants [22]. Partial thermalization describes peripheral RHIC 
data very well and, moreover, allows one to describe pp and AA collisions in the same model 



Just as with ^ we can use the same correlation function (19) to study flow fluctuations 
( |10] ) and therefore examine the effects of event anisotropy. If flow fluctuations arise only from 
changes in the shape of the collision area, then their governing correlation function would predict 
{Spti6pt2) = and have almost no dependence on collision energy. To start we calculate (10) 
for 200 GeV Au+Au collisions. Using this along with the experimental definition 



Vn{2}^ - Vn{iY 



(20) 



and taking fn{4} ~ (vn) we calculate V2{2} as the solid line in Fig. |3j The dashed line is ^ 
calculated from blast wave and fiow fluctuations (10) determine the difference between the two 
curves. Notice that in central collisions the circular event geometry drives (^2) to zero. Here 
flow fluctuations are largest and completely determine f2{2}. 

As with multiplicity fluctuations and pt fluctuations, the correlation strength TZ sets the scale 
for flow fluctuations. In our model, the increase in growth of flow fluctuations (10) is coupled 
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Figure 4. Measured V2{2} (top) 
from STAR [33] and ALICE [351 IS] 

compared to calculations using ^ 
and (10). Same for i'2{4} (bottom) 
computed from ^ in the blast 
wave. 
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Figure 5. Same calculation as in 
Fig. |4]but for the fourth order flow 
coefficients f4{2} (solid) and f4{4} 
(dashed). ALICE data [35]. 



Figure 6. Flow fluctuations in 
terms of the coefficient of variation 



for elliptic flow (21) compared to 
STAR data l33l. 



to the growth in the saturation scale Qj via (18). Comparisons of V2{2} and f2{4} to different 



collision energy measurements are shown in Fig. |4| It is also informative to make comparisons 
to V4 as shown in Fig. [5j The effect of the change in blast wave parameters with collision energy 
is apparent from the separation of the solid lines in the bottom panel of Fig. [4] and the dashed 
lines in Fig. [5] The importance of the Glasma contribution is similarly evident by the change 
is separation between fn{4} and fra{2} with collision energy. 
STAR has also measured f 2 fluctuations in the form of 



(21) 



resembling the so-called "coefficient of variation" defined as the standard deviation divided by 
the mean. Care should be taken here since the definition of ci^, Eq.(20), is not strictly the 



variance. In Fig. 6 we compare our calculation of (21) to measurement. Calculations at RHIC 



energies seem to agree reasonably well and we include the calculation for Pb+Pb 2.76 GeV as 
a prediction. 

In Fig. [7] we show V3 fiuctuations from Glasma. We further emphasize that the energy 
dependence in Fig. [Tjis in good accord with data, supporting the Glasma scaling with Q^. The 
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Figure 7. Triangular flow fluctu- 
ations, £7^3, Eq. ([To]). Data points 
are computed from f3{2} and f3{4} 
measurements [371 EH] • 



Figure 8. The triangular flow co- 
efficient, vs{2} compared to STAR 
and ALICE data |3ZIES]. 
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Figure 9. Near side ridge 
amplitude calculation from Glasma 
source correlations. Experimental 
data is from (STAR) liSl [29l EO] 
and (ALICE) [38]. 



shape with centrality reflects contributions not only from Glasma, but also our parameterizations 
of the average reaction plane eccentricity e and our choice the flux tube distribution p^^, . 



To use (10) and (20) to calculate f3{2}, we must come to grips with the fact that our 
blast wave parametrization assumes f3{4} = 0. While this seems to be the case for the STAR 
measurements, ALICE has measured a non-zero t'3{4} for Pb+Pb collisions at 2.76 TeV. To 
correct for this possible discrepancy, we can provide an ad-hoc parameterization of W3{4} and 
use (10) and (20) to calculate U3{2}. Agreement shown in Fig. [sjis reasonable. As a preferable 
alternative, we compare our calculated an to fluctuations extracted from ALICE and STAR 
measurements in Fig. [Tj 

The characterization of the ridge in terms of Fourier flow coefficients can prove a valuable 
check for models describing flow and its fluctuations with long range correlations. Relative 



contributions of flow and flow fluctuations to the ridge decoupled using Eq.(12). Though it 
would be more informative to extract an directly from ridge measurements, we can compare 
the flow-subtracted near side ridge amplitude to integration of ([2]), suitably normalized to the 
mixed event background, see Refs.[2l 13 16] for details. 

Our calculation of the flow-subtracted near side ridge amplitude can be seen in Fig. [9] 



The dimensionless coefficient k in (18) is set by the change in amplitude of the 200 GeV 



Au-t-Au measurement and remains fixed for all other calculations. Like the other observables, 
the change in collision energy is determined by TZ and therefore Q^. The agreement with 
data is quite reasonable. Commonality of the energy dependence between the ridge and fiow 
fiuctuations should not come as a surprise. Measurements of f2{2} at ALICE using particle pairs 
with and without a rapidity separation only show noticeable difference in peripheral collisions 



[Ml [39]. The remaining difference between V2{2} and V2{4} suggests (20) results from a long 
range correlation. However, long range correlation explanation of the energy dependence of 
multiplicity and momentum ffuctuations does come as a surprise. Jets and resonance decays, 



for example, should have significant effects on these observables. Further investigation into the 
relative contributions of long and short range correlations can indicate if indeed long range 
correlations dominate fluctuation observables at these energies. 



4. Discussion 

In this work we study the connection between correlations and fluctuations of multiplicity, 
momentum, and harmonic flow coefficients. We first argue in Sec. [2] that given the existence 
of genuine two-particle correlations in the hadron pair distribution, Eq. ([T]), multiplicity 
fluctuations, momentum fluctuations, flow fluctuations, and the ridge are all linked by a common 
scale TZ. We then consider contributions to the genuine correlations due to fluctuating or lumpy 
initial parton distributions; correlated partons emerge from the same source/lump. We take 
the origin of initial state parton density lumps from Glasma flux tubes and present our model 
calculations in Sec. [31 

We focus on two important features of the data: the agreement with the change in collision 
energy and the link between the correlation and fluctuation observables. Long range correlations 
such as those measured in the ridge indicate that the correlations originate at early times in the 
collision. The connection between the ridge and flow fluctuations is evidence that the difference 
between Un{4} and i'n{2} is also a long range phenomenon. Given that flow fluctuations result 
from the existence of parton density lumps in the initial state, one can compare the behavior 
of flow fluctuations to that of other consequences of the lumpy initial state such as multiplicity 
and momentum fluctuations. The correlation strength TZ sets a common scale for all of these 
observables and controls much of the centrality and collision energy dependence. Theories of the 
initial state can therefore be constrained with measurements of TZ. In this work we have chosen 
a CGC-Glasma based description of TZ and found reasonable agreement in the centrality and 
collision energy changes in all of the mentioned observables. This suggests that saturation based 
descriptions of the initial state are consistent with at least mid-central to central collisions. 
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